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' We examine the Detrended Fluctuation Analysis (DFA), which is a well-established method for the detection 

, of long-range correlations in time series. We show that deviations from scaling that appear at small time scales 

. become stronger in higher orders of DFA, and suggest a modified DFA method to remove them. The improve- 

ment is necessary especially for short records that are affected by non-stationarities. Furthermore, we describe 
^ , how crossovers in the correlation behavior can be detected reliably and determined quantitatively and show how 

■ several types of trends in the data affect the different orders of DFA. 
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I. INTRODUCTION 

In recent years the Detrended Fluctuation Analysis (DFA) invented by Peng et al. has been established as an important 
. tool for the detection of long-range (auto-) correlations in time series with non-stationarities. It has successfully been applied to 



such diverse fields of interest as DNA [|3]-|], heart rate dynamics [p[-p3|], neuron spiking 1 16,r7|], human gait [|lO| 



[18], long-time 

weather records Jl^-pl]], cloud structure [ !2[^], economical time series [p4[-p7||, and even solid state physics Pq, 29 ] . While the 
spectral analysis (Fourier transform) and the recently developed wavelet transform modulus maxima (WTMM) method [^0[-p2[] 
analyze the time series directly, the DFA is based on random walk theory [B3,p4|, similar to the Hurst rescaled-range analysis 



> 

^ [ 135| ] (see also [|36|]) and similar to another method based on wavelet transform used e. g. in [[19|^0|]. Since the time series are 
summed in the methods based on the random walk theory (including the DFA), the noise level due to imperfect measurements 
in most records is reduced. The fano factor [|7|] and the allan factor [^], see also [H], have been employed in a similar context, 
but these methods do not remove trends in the data. 
, For the reliable detection of long-range correlations, it is essential to distinguish trends from the long-range fluctuations 
^ intrinsic in the data. Trends are caused by external effects - e. g. the greenhouse warming and seasonal variations for temperature 
O records - and they are usually supposed to have a smooth and monotonous or slowly oscillating behavior. Strong trends in the 
I ^ I data can lead to a false detection of long-range correlations if only one (non-detrending) method is used or if the results are not 
carefully interpreted. It is the advantage of the DFA that it can systematically eliminate trends of different order (like the method 



^-H based on wavelet transform that has been applied e. g. in 1 19,20p). This way we can gain insight into the scaling behavior of the 
K*" ■ natural variability as well as into the trends in the considered time series. 

In this paper, we study systematically different orders of the DFA technique, that allow to eliminate different orders of trends. 
The paper is organized as follows: In Section 2 the method is described. In Section 3 we suggest a straightforward extension of 
the DFA that eliminates DFA specific deviations from scaling at small time scales. We describe how crossovers in the observed 
long-range correlation behavior can be detected and detail how the crossover time can be determined reliably. Finally, we show 
T— I how several types of trends in the data affect the different orders of DFA. We summarize the results in the forth section of the 
'"^ paper. 



II. LONG-RANGE CORRELATIONS AND THE DETRENDED FLUCTUATION ANALYSIS (DFA) 



^ ' We consider a record (xi)ofi = 1, . . . ,N equidistant measurements. In most applications, the index i will correspond to the 
Q time of the measurements. We are interested in the correlation of the values Xi and Xi+s for different time lags, i. e. correlations 
^ over different time scales s. In order to get rid of a constant offset in the data, the mean {x) — X^ili usually subtracted, 
Xi = Xi — {x). Quantitatively, correlations between x-values separated by s steps are defined by the (auto-) correlation function 

C{s) = (xiXi+s) = y^XiXi+s- (1) 



X 



=1 



If the Xi are uncorrelated, C(s) is zero for s > 0. Short-range correlations of the Xi are described by C(s) declining exponen- 
tially, C{t) ~ exp(— s/sx) with a decay time Sx- For so-called long-range correlations C{s) declines as a power-law 

C(s) oc s-T (2) 

with an exponent < 7 < 1. A direct calculation of C{s) is usually not appropriate due to noise superimposed on the collected 
data Xi and due to underlying trends of unknown origin. For example, the average {x) might be different for the first and the 
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second half of the record, if the data are strongly long-range correlated. This makes the definition of C(s) problematic. Thus, 
we have to determine the correlation exponent 7 indirectly. 

Often experimental data are affected by non-stationarities. Such trends have to be well distinguished from the intrinsic 
fluctuations of the system in order to find the correct scaling behavior of the fluctuations. This task is not easy, since e. g. 
subtracting some kind of moving average with a certain bin width a would artificially introduce the time scale a into the data, 
thus destr oying a possible scaling over a wider range of time scales. Hurst rescaled-range analysis and other non-detrending 
methods [ p7| , p8| ] work well if the records are long and do not involve trends. But if trends are present in the data, they might 
give wrong results. Very often we do not know the reasons for underlying trends in collected data and - even worse - we do not 
know the scales of the underlying trends. Detrended fluctuation analysis (DFA) is a well-established method for determining the 
scaling behavior of noisy data in the presence of trends without knowing their origin and shape 

The DFA procedure consists of four steps. In the first step, we determine the profile 



k=\ 



(3) 



of the record (xi) of length iV. The subtraction of the mean (a;) is not compulsory, since it would be eUminated by the later 
detrending in the third step anyway. 
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FIG. 1. Illustration of the detrending procedure in the Detrended Fluctuation Analysis (DFA). For two segment lengths (time scales) 
s = 100 (a) and 200 (b), the profiles Y{i) (dashed lines; defined in Eq. (|3|)), least square quadratic fits to the profiles (solid lines), and the 
detrended profiles Ys (i) (dotted lines) are shown versus the index i. 

In the second step we cut the profile Y{i) into Ns = \^ I s\ non-overlapping segments of equal length s (see Fig. Since the 
record length N need not be a multiple of the considered time scale s, a short part at the end of the profile will remain in most 
cases. In order not to disregard this part of the record, the same procedure is repeated starting from the other end of the record. 
Thus, 2Ns segments are obtained altogether 

In the third step we calculate the local trend for each segment by a least-square fit of the data. Then we define the detrended 
time series for segment duration s, denoted by as the difference between the original time series and the fits, 

Ys{i)=Y(i)-Vv(i), (4) 

where is the fitting polynomial in the i/th segment. Figure |l| illustrates this procedure for s = 100 and 200. In the example 
quadratic polynomials are used in the fitting procedure, which is characteristic of quadratic DFA (DFA2). Linear, cubic, or 
higher order polynomials can also be used in the fitting procedure (DFAl, DFA3, and higher order DFA). Since the detrending 
of the time series is done by the subtraction of the fits from the profile, these methods differ in their capability of eliminating 
trends in the data. In ?ith order DFA, trends of order n in the profile and of order ?i — 1 in the original record are eliminated. 
Thus a comparison of th e resu lts for different orders of DFA allows to estimate the strength of the trends in the time series, as 
will be shown in Section [II C. 
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In the fourth step, we calculate - for each of the 2iy5 segments - the variance 

= = 1 ^ Y^[{j. - 1). + ^] (5) 

1=1 

of the detrended time series Ys{i) by averaging over all data points i in the z^th segment. At last we average over all segments 
and take the square root to obtain the DFA fluctuation function 



F{s) 



2Ns 

2/ 
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n 1/2 



(6) 



For different detrending orders n we obtain different fluctuation functions F(s), which we then denote by F^") (s). By construc- 
tion, (s) is only defined for s > 7i + 2. We are interested in the s-dependence of i^*^"^ (s). It is apparent, that the variance 
will increase with increasing duration s of the segments. If the data (x^) are long-range power-law correlated [see Eq. (^], the 
fluctuation functions F'"^ (s) increase by a power-law 

F(")(s)(xs" (7) 

for large s values, where the fluctuation exponent a is related to the correlations exponent 7. A direct derivation of Eq. (Q) for 
large s is given by | ]3^ ] in the appendix. For data without trends and zero offset, Xi — Xi and Ys{i) = Y{i) for i < s. Then the 
mean-square displacement in each segment n can be calculated [with Eqs. ^ and ([l|)]: 

+ ( g x,Xk)=i{x^)+Y. C{\k-j\) = i{x^) + 2Y^{i-k)C{k). (8) 
k=i I \ k^j I k^j k=i 

For large i the second term can be approximated [with Eq. (0)]: 



i — l i „j 2 — 1 

^ C{k) r^^k-'' i k-'<dk - i^^'< and ^ kC{k) 

1 1 I 1 1 I 1 



i^-T. (9) 

k=l fe=l •'-^ fc=l 



If the data are long-range power-law correlated with < 7 < 1, this term will dominate for large i, giving 

(y2(,))^,2-^. (10) 

Thus, the mean-square displacement {Y'^{i)) of the profile increases faster than linearly in i, which corresponds to superdiffu- 
sion. A similar approximation for F^'^\s) using Eq. (|l^ finaUy leads to 

for large time scales s (see the appendix of js^ ] for an exact derivation for DFAl). Thus, comparing Eqs. (Q) and ([ll|), we find 

a = 1-7/2 for 0<7<1. (12) 

If the data are uncorrelated or short-range correlated [C(s) decays exponentially or 7 > 1 in Eq. (^], the first term in Eq. 
will dominate for large i, and we find {Y^{i)) ^ i (corresponding to regular diffusion) and hence ^("^(s) ^ s^/^. Thus, the 
fluctuation exponent a = 1/2 indicates the absence of long-range correlations. 

Practically, we can plotF(")(s) as a function of s on double logarithmic scales to measure a by a linear fit. For uncorrelated or 
short-range correlated data, we expect a = 0.5, while a > 0.5 indicates long-range correlations. In this case we can determine 
the correlation exponent 7 by measuring the fluctuation exponent a. Figure ^ shows two examples for the application of the 
DFA method to long-range correlated data [Fig. ||(a,b)]. An example for uncorrelated artificial data [Fig. ||(c)] is also shown, 
confirming a = 1/2. 
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FIG. 2. Detrended Fluctuation Analysis (DFA) of long-range correlated data with (a,d) 7 = 0.2 (a — 0.9), (b,e) 7 — 0.7 (a = 0.65), 
and uncorrelated data (c,f) with a = 0.5. The scaled common DFA fluctuation functions _F'"'(s)/s^/^ are plotted versus the time scale s in 
(a,b,c), while the corresponding modified DFA fluctuation functions Fj^J^{s) are shown in (d,e,f). The symbols correspond to the different 
detrending orders n, DFAl (□), DFA2 (O), DFA3 (A), DFA4 (y), DFA5 (O), and DFA6 (+). The artificial long-range correlated data have 
been generated by the Fourier transform method, see e. g. [ ^p^ . The deviations from scaling occurring in F(s) for small time scales s in 
parts (a,b,c) are drastically reduced by dividing by the correction function obtained from shuffled data in the parts (d,e,f). The results have 
been obtained by averaging over 100 artificial series of length A'^ — 200, 000 for each part. 



III. RESULTS 



A. The Correction Function and an Modified Version of the DFA 

Figures ^a,b,c) show that small deviations from the scaling law Eq. (0), i. e. deviations from a straight line in the log-log plot, 
occur for small scales s. These deviations are intrinsic to the usual DFA method, since the scaling behavior is only approached 
asymptotically. The deviations limit the capability of DFA to determine the correct correlation behavior in very short records 
and in the regime of small s. DFA6, e. g., is only defined for s > 8, and significant deviations from the scaling law occur 
even up to s « 30. They will lead to an over-estimation of the fluctuation exponent a, if the regime of small s is used in a fitting 
procedure. Previously, an attempt has been made to improve the scaling by modifying the prefactor in the definition of i^(s) in 
Eq. 1^], but it is valid for DFAl only and it still contains some approximations. Here, we suggest a different approach, which 
can also be applied for higher order DFAs. 

For long artificial series of uncorrelated (a ~ 1/2) and long-range correlated data (with a — 0.9) we have determined the 
deviations of F*^"^ (s) from the expected scaling behavior Eq. (^. Figure ^ shows the results for the correction function 

where, again, n denotes the DFA detrending order, and (. . .) denotes the average over different configurations. Practically, s' 
has to be large (s' > 50), but it must remain significantly smaller than the record length N; s' « A^/20 seems to be a reasonable 
number. If we divide the DFA fluctuation functions F{s) by the corresponding correction function A'!"-* (s), the deviations from 
scaling for small s are eliminated. The crucial point is that the correction function ifi"^ (s) depends only weakly on a (see 
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Fig. |3]). Therefore, practically, the correction function for uncorrected data, k'^^}^{s), can be used in all cases. k'^^^^{s) can be 
obtained most easily by analyzing the corresponding shuffled data, where all long-range correlations have been destroyed by 
randomly shuffling the record of measurements. 
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FIG. 3. The DFA correction functions K^\s) for uncorrelated data (a = 1/2, small filled symbols) and for long-range correlated data 
(q = 0.9, large open symbols) are shown versus s. The symbols correspond to the different detrending orders of DFA (see Fig. |^ for the 



definition) and the functions are shifted by multiples of 1.3. The dashed lines indicate the asymptotic behavior. K^\s) 1. In the 



numerical procedure, (s) has been obtained by averaging over 100 artificial series of length A'^ = 200, 000. The long-range correlated 
series have been generated by the method of Fourier transform, see e. g. 

Thus, in order to improve the scaling of the DFA fluctuations on short scales s, we suggest to divide F*^"^ (s) by iirj"2('S) 
obtain the modified fluctuation function 



, ^ F(")('s^ , > (fF^^Lfs'll^^l/Z .,1/2 



(14) 



according to Eq. (p3|). Here, ([^^Jhuff (s)]^)^^^ denotes the usual DFA fluctuation function [defined in Eq. (^] averaged over 
several configurations of shuffled data taken from the original record [x-i) under consideration, and s' k, N/20, again. 

In Fig. ^ the numerical results for Fj^^^{s) are compared to those for the common F^"-^ (s). The improvement of the scaling 
behavior for small s can be seen clearly by comparing the right and left parts of the figure. Only very weak deviations from the 
expected power-law dependence remain for strong long-range coiTelations [for large a, see Fig. §(d)]. The suggested modified 
method does only need additional computation time, but the programming effort is not significantly higher. Thus, it can easily be 
used in most applications. The improvement is very useful especially for short records or records that have to be split into shorter 
parts to eliminate problematic nonstationarities, since the small s regime can be included in the fitting range for the fluctuation 
exponent a. 

Note that there is another advantage in dividing by the DFA fluctuation function for shuffled data: If the distribution of the 
x- values in the record is very broad, e. g. similar to a Levy distribution, and the second moment (a;^) diverges, systematic devi- 
ations from the expected scaling behavior will occur already for uncoiTelated data (see e. g. [0]). Since these deviation are not 
eliminated by shuffling the data, they will cancel out in Fj^^^{s). Thus, the modified DFA method indicates the correct correla- 
tion behavior also in presence of broadly distributed data, where the common DFA fails to distinguish long-range correlations 
from deviations caused by broad distributions. 



B. Determination of Crossovers 



Frequently, the correlations of recorded data do not follow the same scaling law for all time scales s, but one or sometimes 
even more crossovers between different scaling regimes are observed. For example, the data might become uncorrelated on large 
time scales s > s^. In such a case, it would be useful to extract the crossover time from the data also by means of DFA. In 
order to do this, we have to investigate how crossovers in the correlation properties show up in the DFA fluctuation functions 
with different orders of detrending. 
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Artificial time series with a well-defined crossover at Sx are most easily generated in a modified Fourier transform procedure: 
The power spectrum P{f) of an uncorrected random series is multiplied by {f / fx)^^ with f3 — 2a — 1 for frequencies 
f > fx ^ ^/sx only. The series obtained by inverse Fourier transform of this modified power spectrum exhibits power-law 
correlations on time scales s < Sx only, while the behavior becomes uncorrected on larger time scales s > Sy. An inverse 
crossover with long-range correlations only for s > Sx and uncorrected behavior below Sx is obtained in a similar way, if we 
multiply the power spectrum by (///x ) for low frequencies f < fx only. Note that there is an alternative way to generate 
series with a crossover in the correlation behavior: We can divide the original long-range correlated series into segments of 
length Sx and shuffle the segments. This way, all correlations for s > Sx are destroyed, but the correlations within the segments, 
i. e. for s < Sx, are preserved. With this method the crossover regime turns out to be much broader, so we use the Fourier 
transfrom method in this paper 
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FIG. 4. Modified DFA of artificial series with a crossover from long-range correlated behavior [(a): 7 — 0.2, a — 0.9; (b): 7 = 0.7, 
a = 0.65] for s < Sx = 200 to uncorrelated behavior for s > Sx- The scaled modified DFA fluctuation functions F^J^{s)/ s^''^ are 
plotted versus the time scale s for DFAl to DFA5. In part (c) the results for data with an inverse crossover from uncorrelated behavior (for 
s < Sx ~ 200) to long-range correlated behavior (7 — 0.2, a — 0.9 for s > Sx) are shown. For comparison, part (d) shows the usual 
DFA fluctuation functions F{s) for the artificial series already considered in part (b). The long-range correlated series with crossover have 
been generated by the modified Fourier transform method described in the text, and the results for 100 series of length A'^ = 200, 000 have 

been averaged. The symbols correspond to the different orders of DFA (see Fig. □ for the definition). The dotted fits in part (a) illustrate 

(3) 

the procedure used to determine the observed position of the crossover s ^ , while the dashed vertical lines indicate the real position of the 
crossover Sx ~ 200. 



Figure ^ shows the results for the modified DFA fluctuation function i^i„" (j(s), see Eqs. (^ and (p4[), for artificial data with a 

crossover in the correlation behavior. The crossover is clearly visible in the results, but it occurs at times Sx""* that depend on the 
detrending order n and that are different from the original s x used for the generation of the data. This systematic deviation is 
most significant in the DFAn with higher order detrending. It occurs independent of the values of the fluctuation exponents and 
independent of the direction of the crossover (from small to large exponents or vice versa). 

The deviation of the crossovers is systematically investigated in Fig. pi where the position of the original crossover Sx used 

(n) 

for the data generation is plotted versus the position of observed crossover s\, for DFAl to DFA5. The plot can be used to 

(n) 

determine the real crossover position Sx from the s\, estimated with the modified DFAn. If several orders of DFA are used in 
this procedure, several estimates will be obtained which can be checked for consistency or used for an error approximation. 
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FIG. 5. Systematic investigation of tiie crossover detection with tlie modified DFAl to DFA5. For artificial data witii a crossover from 

In) 

long-range correlated behavior (7 — 0.4, a = 0.8) for s < Sx to uncorrelated behavior for s > Sx the crossover times have been 
determined from the intersection of linear fits done on both sides of the crossovers [as illustrated in Fig. ^a)]. The original crossover time 
Sx is varied and plotted versus s'^^' for all five types of DFA. Note that the results also hold for inverse crossovers and for other values of a, 
although the fitting procedure will be less accurate if a is close to 0.5. They also hold for the usual DFA for Sx > 200, since K'^'J^{s) ~ 1 

for large s. The dashed line corresponds to s5^"' = Sx, and its position shows that the estimated crossovers are always larger than the real 
Sx . The long-range correlated series with crossover have been generated by the modified Fourier transform method described in the text. For 
each point the results for 200 time series of length A'^ — 100, 000 have been averaged, and the error bars have approximately the size of the 
symbols. The symbols correspond to the different orders of DFA (see Fig. for the definition). 
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FIG. 6. Investigation of the trend elimination and detection capability of the conventional Fluctuation Analysis (FA, symbol — ) and the (not 
modified) DFAl to DFA5 for trends Ax^ with integer powers p. For artificial long-range correlated data (7 = 0.7, a — 0.65) with added (a) 
linear, (b) quadratic, (c) cubic, and (d) 4th order trends, the DFA fluctuation functions F'^"\s) are plotted versus the time scale s (A = lO**, 
p — 1,2,3, 4). The trends are completely eliminated if n > p. This feature allows to determine the (integer) order p of the trends using DFA. 
For n < p the trends lead to an apparent crossover at high s values as discussed in the text. Only one series of length A'^ = 100, 000 has been 
considered for each part of the figure. The symbols correspond to the different orders of DFA (see Fig. pifor the definition). 



C. Monotonous Trends 



Records from real measurements are often affected by trends, which have to be well distinguished from the intrinsic fluctu- 
ations of the system as discussed in Section 2. To investigate the effect of trends on the DFA fluctuation functions, we have 
generated artificial series (x^) with smooth monotonous trends by adding polynomials of different power p to the original record 
{xi) generated with the Fourier transform method: 



AxP 



with X — i/N. 



(15) 



Figure ^ shows the effect of trends of different (integer) order p on the FA and DFA fluctuation functions F'^")(s). For the 
conventional Fluctuation Analysis (FA), the variance F?^{i^) for segment v of scale length s is determined by the mean square 
deviation of the (non-detrended) profile Y{i), i. e. Eq. (pb is replaced by 



F^{iy) = [Y{vs) - Y{vs -s + l)f 



(16) 
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FIG. 7. Investigation of the trend elimination and detection capability of the modified DFAl to DFA5 for trends Ax^ with non-integer 
powers p. For artificial long-range correlated data (7 = 0.7, a = 0.65) with added trends (as listed in the figure), the scaled modified DFA 
fluctuation functions Fj^^^{s) /s^^^ are plotted versus the time scale s. As expected, the trends are best eliminated by the higher order DFAs, 
but an apparent crossover is still observed for large time scales s. The position of the crossover and the slopes of the curves above the crossover 
can be used to determine the strength A and the exponent p of the trends (see text). The results for 100 series of length TV — 200, 000 have 
been averaged. The symbols are the same as in Fig. 
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This procedure leads to similar results as the DFA, but trends are not eliminated. In Fig. ||the strong trends completely conceal 
the long-range correlations of the time series. The slope atrcnd = 1 is observed, which is the maximum for the FA, characteristic 
of strong trends (or non-stationary time series). For the DFA, the trends in the data can lead to an artificial crossover in the scaling 
behavior of F'^")(s), i. e. the slope a is increased for large time scales s. The position of this artificial crossover depends on the 
strength A and the power p of the trend. Evidently, no artificial crossover is observed, if the detrending order n is larger than 
•p. Thus, the order p of the trends in the data can be determined easily by applying the different DFAn. If -p is larger than n, an 
artificial crossover is observed, and the slope ationd in the large s regime strongly depends on n. Thus, the artificial crossover 
can be clearly distinguished from real crossovers in the correlation behavior, which would result in identical slopes a and rather 
similar crossover positions for all detrending orders n (see Fig. 

Figure ^ shows the same type of results as Figs. ^, but for the modified DFA and for trends Ax^ with non-integer power p. 
Here, strong trends are not completely eliminated by the DFA and can dominate the behavior of the fluctuation functions on very 
large time scales s even if n > p. If the trends are not too weak [as for the DFA5 in Fig. ^a)] or too strong [as for the DFAl in 
Fig. ^c,e,f)] an artificial crossover occurs. The positions of the artificial crossover and the slopes Q!|"cnd ^'^ '^^^ large s regime 
can be used to determine the trend parameters A and p. Two main rules can be deduced from the results shown in Fig. ^ and we 
have confirmed these rules also for other parameters a, A and p. (i) If a trend-related crossover is observed and the detrending 
order n is sufficiently small (n < p + 0.5), the position of the artificial crossover depends only on A, but not on p or a [compare 
Figs. ^c,e,f)]. This permits to determine the strength A of the trend for real data by comparison to the results for artificial data 
with known trend strength [see Figs.|7|(a,b,c)]. (ii)If a trend-related crossover is observed, the slope ck|"jj^j for large s (above 
this crossover) is the minimum of n + 1 and p + 1.5 [see Fig. ^c,d,e)]. This permits to determine the order p of the trend, if n 
is chosen sufficiently large. The combination of both rules allows to determine A and -p if the trend is strong enough and several 
appropriate DFA orders n are employed. Note that the high order DFAs tend to become numerically unstable on small scales s 
for very strong trends and weak fluctuations [see DFA7 for s < 30 in Fig. 0(c,e,f)]. 




FIG. 8. Modified DFA of long-range correlated random data (a = 0.65, unit variance) with additional oscillatory trends y4sin(27ri/) of 
different frequencies /, (a) / = 0.1 and (b) / = 0.001, and intensities A ranging from A = 0.02 to ^4 = 2.56. The modified DFA fluctuation 
functions f£^Q(j(s) are plotted versus the time scale s. For other orders of the DFA we get similar results. Only for strong oscillatory trends 
the expected scaling behavior F^^^(^s) ~ s^'^^ is disturbed. This shows that DFA is quite robust against oscillatory trends, especially against 
those with rather high frequencies. The results for 100 series of length A'^ = 200, 000 have been averaged for each curve. 



D. Oscillatory Trends 

Sometimes the trends superimposed on the records of real measurements are not monotonous. A quite common shape of trends 
is slow oscillatory behavior. It can be either rather periodic and regular, e. g. for the seasonal trend in temperature fluctuations, 
or rather irregular, e. g. for trends due to the southern oscillation (El Nino) effect. In both cases, the effects of the trends have 
to be distinguished from the intrinsic fluctuations of the system under consideration, and strong oscillatory trends can lead to a 
false determination of the correlation behavior. In particular, it is of great interest to know, how strong the oscillatory trends may 
become until they start to disturb the correlation analysis. 

In order to investigate the effect of oscillatory trends on the results of the DFA, we have constructed artificial records of 
long-range correlated random numbers with superimposed oscillatory trends of different frequency and strength. The results for 
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the modified DFA fluctuation function F^^^{s) are shown in Fig. ^. The figure shows that the DFA is more sensitive to slowly 
varying trends, while quickly oscillating trends disturb the scaling behavior of the results much less. Thus, for example, the 
seasonal trends in temperature records have to be removed prior to the analysis (e. g. by subtracting a daily mean temperature), 
while the modulating effect of the breathing does not significantly disturb the analysis of heartbeat records. 
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